INVARIANCE OF THE DISTRIBUTIONAL CURVATURE 
OF THE CONE UNDER SMOOTH DIFFEOMORPHISMS 



J.A.VlCKERS AND J.P.WlLSON 

- - - Abstract. An explicit calculation is carried out to show that the distributional curvature of a 2-cone, calcu- 

lated by Clarke et al. (1996), using Colombeau's new generalised functions is invariant under non-linear 
coordinate transformations. 

0^ 



X 



1. Introduction 



Recently, Colombeau's theory of new generalised functions (Colombeau, 1984) has been applied to 
give a distributional interpretation to physical quantities in General Relativity that are calculated by 
non- linear processes, such as the distributional curvature of cosmic strings (Clarke et al, 1996), the 
Ultra-relativistic Riessner Nordstr0m field (Steinbauer, 1997) and the energy-momentum tensor of the 
^ ■ Kerr solution (Balasin, 1997). One of the major problems of using the full Colombeau algebra, which 

OO I undermines its usefulness in a covariant physical theory such as Ceneral Relativity is that it is not of a 

' coordinate invariant construction; that is given a diffeomorphism /i : O — > f2' between two open sets, there 

is no natural way of defining a map //* : £(Q') — > which preserves the concepts of moderate and 

' null functions, and commutes with the canonical embedding (smoothing convolution) l : V^fl) £{Q). 

00 ! In this paper we shall return to the calculations of Clarke et al. (1996) and we shall explicitly show 

0^ ' that the distributional curvature is independent of any non-linear C°° coordinate transformation that 

^ I we might carry out on the metric. This will provide evidence for a covariant formalism of Generalised 

O^' functions 
U 

2. The distributional curvature of cosmic strings 



The distributional curvature of the cone with a deficit angle of 27r(l — A), whose metric may be written 

(1) 



as 

gab = i(l + A'^)5ab + i(l - A^)mab 



2 > 



\ x'^+y'^ x^+y^ / 

was recently calculated by Clarke et al. (1996) by using Colombeau's theory of generalised functions 
to overcome the problem of assigning distributional interpretations to products of distributions. It was 
shown to be 

i?V5 = 47r(l-A)<5(2)(a,,y) 

This result implies that a thin cosmic string, whose exterior metric is 

ds^ = -dt^ + dr^ + A^r^d(j}'^ + dz^ 

has a mass per unit length of 27r(l — A)/ A. 

The process used to calculate the curvature first involves an embedding of the metric gab into £m{^^)] 

L ■ gab>-^ g^b 

g7b{^,x,y)^ / gabix + e^,y + er])^{^,'q)d^dr] 
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The generalised function Ricci curvature R^/g is then is calculated from the metric g^i, in the usual 
manner, and finally it is shown that this curvature is equivalent to the distribution 47r(l — A)5^^^ in the 
sense of weak equivalence; that is that for each ^ e '"'"w^'' 



lim / Rj^^dxdy = "^{0,0) 



for $ G Aq{M?), for large enough q e N. 

The theory of General Relativity is however a covariant theory, so we would like these results to make 
sense in the context of coordinate invariance; that is given a C°° diffeomorphism; 

{x',y') =fJ,{x,y) 

the following diagram commutes; 

9U Zb ^ R'VW 47r(l-A)^(2) 

gab gZb > RVg 47r(i - 

where /x* denotes the appropriate tensor or density transformation laws 

Since the full Colombeau algebra is not a coordinate invariant construction and the embedding t docs 
not commute with /z, it will only make sense to talk about coordinate invariance at the level of distribution 
theory; this will mean that we must follow the approach of carrying out the coordinate transformation 
on the metric g^^ to form a new metric gab, explicitly calculate the distributional Ricci scalar density 
for the new metric and show that it is transforms as a scalar density of weight +1 from the Ricci scalar 
density that corresponded to g'^f^. 

Wc first apply C°° diffeomorphism /i : ^ to the primed version of the metric (1); The fact that 
Colombeau's theory is manifestly invariant under linear coordinate transformations means that without 
loss of generality, it may be assumed that /U is of the form 

x' = x + f{x,y) 
y' ^y + g{x,y) 

where / and g are C°° and 0{x^ + y'^)- 

On applying (2) to (1) we arc able to write 

gab{x, y) = i(l + A2)Za6(x, y) + ^{1 - A'^)mab{x, y) 

where 

_ dx^dx^ 
~ dx- dx^ 



and 



''ab 



{x+!f-{y+gY 2(x+/)(y+3) 
2(a;+/)(y+g) {x+!?-{v+gf 



V 



Next we shall smooth gafc with a kernel $ e ^o(IR^) having a support in a radius of 

i?o = {(a;2 + 2/2)V2 | \^[x,y)\>Q]. 
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The functions lab, consisting of products of the smooth Jacobian components, will be C°° and so may be 
identified with their smoothings in £m(R^)- 

The matrix niab, containing the singular contributions to the metric, may be expressed as 



rriab = 1 



,xf + yg , P + g' 



2xy+f2l 

3,2_|_y2 



'2xy+fi2 



where the fab{x,y) are C°° and 0{r^) functions. It may be observed that the components of the above 
matrix may be expressed as linear combinations of real and imaginary parts of complex- valued functions 
having the form 

, . ^ e^'-^ + mi {x, y)e-'^^^ + ma (a;, y) + msjx, t/)e^''^ 
1 + mi{x,y)e-'^'-'t' + m5(a;,y) + m6{x,y)e'^^'l' 

where we define the polar coordinates (r, <p) by 



X = r cos ( 



y = r sm ^ 



and the mk{x,y) are C°° and 0{r). Therefore we may calculate niabe by smoothing h{x,y). 
The expression for m{x,y) may be expanded as 



m 



(x, y) = e'^f' + r ^ aue^'^ + rM^, y) 

fc=±l,±3 



(3) 



where the ak are constants and p{x,y) is 0(r°), continuous everywhere and C°° everywhere except at 
r = 0. We may smooth this expansion termwise and obtain estimates for the smoothings. Certainly the 
re^"^ are C°° and may be identified with their smoothings. For other terms, the results in the appendix 
will give 

Ci + C2xle + Czy/e + 0{r^/e^) ii r < eRo 



^2i<p — 



ii + o 



if r > eRo 



' £(C±4 + Cisx/e + C±6y/£ + 0(rV£')) if r < eRo 



re 



±3i. 



1 + 



if r > eRo 



where the Ck are constants. 

The only smoothing that we have not yet had to estimate is that of the remainder term r^p{x,y). 
This term will provide a higher order contribution to the metric and therefore will not require a such a 
delicate estimate in order to be able to differentiate it when we come to calculate the curvature. We may 
write the smoothing of this term as 



^ /.OO />27r 

r'^p{x, y) = / (r^ + 2epcos(V' — (/))+ e^p^) 
Jo Jo 



For r < eRo, we have 

r'^p{x,y) 



dx 



|r2p(x,y)| 



r=0 



r=0 



X p{r COS ^ + ep cos V', r sin ^ + ep sin tp) 
X ^{pcosil), p sin il))p dp dip 



■■e^ [ {X^ + p,'^)p{eX,ep,)^X,p)dXdp = 0{e'^), 
■ 2e / Xp{eX,ep)^{X, p) dXdp 
+ / (X"^ +p'^)p^^{eX,ep)^{X,p)dXdp = 0{e), 
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and similarly 



Therefore 



dy 



|r2p(a:,2/)| 



= Oie). 



r=0 



r'^p{x, y) = M{e) + 0{er) 



where M{e) is an O(e^) function of e alone. 
If however, r > sRq, then 



r^'pi^y) = (^1 + 2^ cos(^ - + 

X p(^r{cos<p + ep/r cos ij)), r [smcj) + ep/r sin'^)) 
X ^{pcosil), p sin ip)p dp dip 

and so by applying the Mean Value Theorem to f{s) = p{x + sX, y + sfj) 

p{x + \,y + n)= p{x, y) + \p,x{x + i\,y + ^n) + HP,y{x + ^X,y + 

where ^ e (0, 1); it may be concluded that 

r'^p{x,y) = r^p{x,y){l + 0{s/r)). 

We now may add up the smoothings in the summation (3) to obtain 

C\ie) + C2(e)f + C3(e)f + ©(rVe^) + 0(re), if r < eRo, 

- 0{rs), if r > e-Ro; 



express the smoothed metric gab as 



9abe = < 



C'ab (e) + Ctl^ ie)x/e + C]^,' {e)y/e + Oir^e^) + 0{rs) if r < sRo 
gab + o(s^)+ 0{re) if r > eRo 



where C^'J (e) are 0(1) functions of s alone, and use it to calculate the Ricci scalar density as 



ReVds 



( 0{l/s'^) + 0{e/r) ifr<sRo 
\ 0{s/r^) if r > eRo 



To obtain a distributional interpretation to this curvature, we shall follow the method of Clarke et al. 
(1996). Much of the analysis is applicable, but this time we shall be working with more crude estimates. 
Given * G ^(R^), we define 

K = supp ^' 
J?K = sup{(a;2+ 2/2)1/2 | y)| > q } 

B, = sup{ {x, y) e I (x^ + t/2)i/2 < I 
By the mean value theorem, we may write for some ^ G [0, 1], that 

/ Re\^{x,y)^{x,y)dxdy = Ii+l2 
Jk 
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where 



h = [ Re ^Mdxdy^ {0,0) 

JK 



dx dy. 



Now 



/2KM1 / Re^/fe rdxdy + Mi ReVFe 



rdxdy 



^ 2nM2 + 2TTM3e^Ro^ + 2TTM4e{\ogRk - log(ei?o)) 

where M„ are positive constants, giving that I2 = 0(£loge). 

To calculate h we let £> = { {x, y) \ (a;^ + y^)^/^ ^ -Rn } C if. and write 

h = I Re^/fsdxdy'^{0,0)+ [ Res/fedxdy^ {0,0) 
Jk-d jd 

The first integral will be 0{e/R]j), converging to zero as e — > 0. For the second integral the Gauss-Bonnet 
theorem is applicable. It is also the case that 



Therefore 



/ Kg^ ds = 2-kA + 0{s/Rd) 

JdD 



h = 47r(l - ^)*(0, 0) + 0{s/Rd). 



It may therefore be concluded that 



£^0 



lim^ / Re^ge{x,y)^{x,y) dx dy = At:{1 - A)^{0,0) 



K 



giving us 



[R^/^]^47r{l-A)6^^\x,y). 



This shows that the distributional curvature associated to the curvature of the generalised metric 
transforms as a scalar density of weight +1. 



3. Conclusion 

We have now shown that the calculations of Clarke et al. (1996) are independent under C°° diffeo- 

morphisms, which suggests that Colombeau's theory could be reformulated in a covariant manner. 

There have been many attempts at amending Colombeau's theory in order to achieve coordinate 
invariance in this sense. One such approach has been to use the simplified algebra (Biagioni, 1990; 
Colombeau, 1992) in which one defines the space base £{n) as consisting of functionals dependent on 
the regularisation parameter e € (0, 1], rather than the smoothing kernel $ G Ao{^), thus avoiding the 
problem of using the spaces AqW^) in the definitions of £m{^) and Af{fl); however this still does not 
make the embedding, which involves $, commute with /i. 

A completely different approach was pursued by Colombeau and Meril (1994) in which they substan- 
tially redefined the kernel spaces Aq{M.'") with weaker moment conditions imposed which does result in 
the ability to construct Q{Cl) in a coordinate invariant manner together with the embedding commuting 
with 12. This formalism enables scalar generalised functions to be constructed, whose transformation laws 
coincide with those of scalar distributions. 

Work is currently in progress (Vickers and Wilson, 1998) into extending the formalism of Colombeau 
and Meril to enable generalised functions to be constructed as multi-index tensors whose transformation 
laws coincide with those of tensor distributions. 
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Appendix. Smoothing of polar functions 

In this appendix we shall smooth the non-regular complex functions e^**^ and re^*"^ and hence obtain 
estimates for the regions r < eRo and r > eRq. For the fine details of calculation the reader is referred 
to Clarke et al. (1996). We begin with the function 

f{x,y) = r'^-ie'(»+i)«* 

where a = 1,2. Its smoothing, with a kernel $ G Aq{M.'^) may be expressed in polar coordinates as 



/($,rcos (p,rsm cj)) 



(re"* + e/9e*'^)"<i>(pcos'(/',/9sin^) 



+ e2p2 _|_ 2epcos{'ip — tj)) 
We may expand $ as a sum of circular harmonics; 

$(r cos (j), r sin </>) = ^ e^"''^ 



pdp 



nez 



so that we may write 



where 



/($, r cos (t),rsm(f> = ^ F„ (r-)e'("+»+i)'^ 
nez 



On integrating out ip we obtain 



+ ep/re 



-itjj 



if n > 



^ -27rr«-i/,7,(-f )" (l- 2V)"$n(p)prfp if n< -(a+1) 



We are interested in estimates for both r < eRo and r > eRo] 
To estimate / for r < sRq we begin by estimating the integrals 



an.feW = 2n p'^+'p'^+'^nip) dp k^O 
Jo 

f-OO 

Mr) = 27r / p'=+V'+'*n(p) dp A: < 1 

>/r/£ 



It is easily seen that 



about r/e = 0; 









l/3n,fc(r) 


-2 


are more 


/3n,- 


2(r) = 27r 


/3n,- 


i(r) = 27r 


/3n 


o(r) = 27r 


/3n 


i(r) = 27r 



rfe+2 
r.fe+2 
rfe+2 



fc > 



fc ^ 3 



/•OO 

/ p$„(p)dp + 0(rVe2) 
Jo 

/■OO 

/ .2^ /.\^._Lr)/'^3/^3N 
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We may use the estimates for a„,fe and (3n,k to obtain estimates for F„ 
For a = 1 we have 

(-l)"fTr {an,n - ^an,n+2^ for n ^ 



for n ^ —2 
for n = — 1 



giving that 



Fn = < 



' 27r/(,°°/9$_2(/9)(i/9 + 0(rVe2) for n = -2 

2n^J^^_i{p)dp + 0{r^/e^) forn = -l 

-27rf /o°°$_3(p)rfp + 0(rV£2) for n = -3 

^ 0{r^/e^) for n > or n < -4 



and for a = 2 we have 



Fn= < 



(-l)";^-^ [an^n - 2f2-a„,„+2 + fi:a„,„+4j for n ^ 

-i-ir^ [f3n,n - 2fJ/3„,„+2 + ^Pn,n+4) for n ^ -3 

£ (2a_i,i - f^a-i.s) + ^^-1,-1 for n = -l 



for n 



giving that 



' s {2tt j^p^ $_3(p) rfp + 0(rV£2)) for „ = _3 
e (27rf /o°°p$o(p) dp + O(rVe')) for n = -2 
£ (-27rf /„°°p$_4(p) dp + Cl(rVe^)) for n = -4 
_ £0(r^/£^) for n ^ — 1 or n ^ — 5 

Therefore we have 

re^'^ f°° 

e2»</' = 27r/ p<i>-2dp + 2TT / 

Jo ^ Jo 

rp-i't> f°° 

-27r / $_3(p)rfp + 0(rV£2) 

^ io 

j,g3j0 = 27r£ / ^-3{p) dp + 27rre*'^ / p $o(p) 

/■oo 

-27rre-''^ / p$_4(p) rfp + O(rVe) 

JO 

On the other hand if r > sRo then 

2'rr«-i /o°° (-^)" (l - ^y^n{p)pdp, if n > 0, 

EL-„ r (-f )" (-^)' <J>n(p)pdp, if -« < n ^ 
0, otherwise. 

On writing the moment conditions, for $ € Aq{M.'^), in polar coordinates as 



Fnir) 



27r / <^o{p)pdp = 1 

27r / p''+i$„(p) dp = 7^ |n| < c < c + n even; 
Jo 
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it is easily seen that 

Fo{r) = l + 0(£«+Vr«+i) 
F„(r) = 0(£«+Vr«+i), n^O 

and hence that 

= e-^ + O 
= re^'"^ + O 
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